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During  the  three  years  of  this  grant,  nine  papers  were  published  or  accepted  for 
publication  in  peer-reviewed  journals.  They  covered  the  following  topics:  (1) 
potential  vorticity  dynamics,  (2)  baroclinic  instability  in  the  presence  of 
forcing,  (3)  baroclinic  wave  structure  and  transitions,  (4)  the  role  of  mountains 
in  .the  evolution  of  mid-latitude  synoptic  disturbances,  (5)  complex  principle 
component  analysis,  and  (6)  numerical  techniques  verses  nonlinear  analytic  methods. 
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I.  INTRODUCTION 


During  the  three  years  of  this  grant,  9  papers  have  been  published  or  accepted  for  publication 
in  quality  scientific  journals  with  high  standards  of  review,  one  paper  has  been  submitted  and 
one  more  will  be  submitted  for  publication  based  on  full  or  partial  support  of  AFOSR  90-0009, 
including  the  work  of  two  Ph.D.  candidates.  The  following  list  also  includes  a  paper  presented  at 
the  AMS  Conference  and  a  report  (partially  supported  by  this  grant)  on  a  software  package  for 
solving  a  general  second  order  parameter  involved  system  of  elliptic  partial  differential  equations 
on  a  two  dimensional  domain  with  Dirichlet  and/or  Neuman  boundary  conditions.  Other  related 
research  under  AFOSR  89-0462  was  reported  in  the  University  Research  Initiative  Final  Report 
supported  by  that  grant. 

Low  Frequency  Oscillations  in  a  Rotating  Annulus  with  Topography,  by  P.  Bernardet,  A.  Butet, 
M.  Dequ4,  M.  Ghil  and  R.  L.  Pfeffer;  J.  Atmos.  Sci.,  47,  3023-3043, 1990. 

Implementation  of  PLTMG  Version  6.0  on  SCRI4D  VAX  with  IRIS  4000  Graphics,  by  Jian  Yu, 
FSU  publication  no.,  FSU-SCRI-90T-124,  August  27,  1990. 

Effects  of  Wave- Wave  and  Wave-Mean  Flow  Interactions  on  the  Evolution  of  a  Baroclinic  Wave, 
by  A.  Barcilon  and  T.  Nathan;  Geophys.  Astrophys.  Fluid  Dyn.,  56,  59-79,  1991. 

Asymmetric  Ekman  Dissipation,  Sloping  Boundaries  and  Linear  Baroclinic  Instability,  by  H.  Weng 
and  A.  Barcilon;  Geophys.  Astrophys.  Fluid  Dyn.,  59,  1-24, 1991. 

Reflection  of  Hydrostatic  Mountain  Waves  from  Spatially  Nonuniform  layers,  by  W.  Blumen  and 
A.  Barcilon;  Geophys.  Astrophys.  Fluid  Dyn.,  58,  25-43,  1991. 

EXSHALL:  A  Turkel-Zwas  Explicit  Large  Time-Step  FORTRAN  Program  for  Solving  the  Shallow- 
Water  Equations  in  Spherical  Coordinates,  by  I.  M.  Navon  and  Jian  Yu,  Computers  in  Geo¬ 
science,  17,  1311-1343, 1991. 

Potential  Vorticity  Index  Vacillation  in  the  1978/79  Winter:  Its  Relation  to  Teleconnection  Pat¬ 
terns,  by  H.-Y.  Weng;  Q.  J.  R.  Met.  Soc.,  118,  327-350, 1992. 

A  Comparison  of  the  Impact  of  Two  Time- Differencing  Schemes  on  the  NASA/GLAS  Climate 
Model,  by  R.  L.  Pfeffer,  I.  M.  Navon  and  X.  Zou;  Mon.  Wea.  Rev.,  120,  1381-1393, 1992. 
The  Effects  of  Viscous  Asymmetry  on  Baroclinic  Equilibrated  States,  by  H.-Y.  Weng  and  A.  Bar¬ 
cilon;  Proceedings  of  the  Eighth  Conference  on  Atmospheric  and  Oceanic  Waves  and  Stability, 
October  1991. 

Nonlinear  Prediction,  Chaos  and  Noise  by  J.  B.  Eisner  and  A.  A.  Tsonis,  Bulletin  of  the  Amer. 
Meteor.  Soc.,  78,  49-59, 1992. 

The  Role  of  Potential  Vorticity  Forcing  in  Topographically  Induced  Instabilities  by  T.  Nathan; 

Trends  in  Atmospheric  Sciences,  edited  by  E.  Menan,  in  press,  1992. 

Low  Frequency  Oscillations  of  Forced  Barotropic  Flow  by  T.  Nathan  and  A.  Barcilon;  submitted 
to  J.  Atmos.  Sci.,  1992. 

Nonlinear  Evolution  of  Baroclinic  Waves  in  the  Presence  of  Zonal  Mean  Periodic  Forcing  by  T. 
Nathan;  to  be  submitted  to  Tellus,  1992. 

Section  II  contains  a  summary  of  accomplishments  under  this  grant  and  section  III  contains 
abstracts  of  the  above  listed  papers. 
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II.  Summary  of  Accomplishments 


1.  PV  Dynamics 

The  goal  of  this  research  was  to  explore  the  conjecture  that  an  index  based  on  the  zonal 
mean  potential  vorticity  (PV)  gradient  on  an  appropriate  isentropic  surface  may  exhibit  a  better 
relationship  with  the  large-scale  circulation,  and  therefore  with  some  global  teleconnections,  than 
traditional  zonal  indices  based  on  zonal  wind  or  surface  pressure.  This  conjecture  was  based  on 
the  knowledge  that  large-scale  atmospheric  processes  are  governed  entirely  by  the  conservation 
of  potential  vorticity  and  that  all  dynamical  Helds,  including  the  vertical  motion,  can  be  derived 
from  the  potential  vorticity  field  on  isentropic  surfaces.  With  this  in  mind,  Dr.  H.  Y.  Weng 
(Quart.  J.  Roy.  Meteor.  Soc.,  118,  327-350,  1992)  defined  a  PV  index  as  the  zonal  mean  mid¬ 
latitude  PV  gradient  on  the  300  K  isentropic  surface  in  the  Northern  Hemisphere  and  used  ECMWF 
FGGE  Ill-b  data  to  investigate  the  evolution  of  that  index  and  its  relationship  with  teleconnection 
patterns  of  the  500-mb  geopotential  height  anomaly.  This  index  fluctuated  with  dominant  periods 
of  the  order  of  two  weeks  and  30-50  days.  It  was  found  that  the  formation  of  blocking  highs  in 
Scandinavia,  Greenland  and  the  Pacific  occurred  at  rather  specific  phases  of  the  evolution  of  the 
PV  index,  especially  as  the  flow  transited  from  low  to  high  values  of  this  index.  The  blocking  highs 
proceeded  in  an  orderly  fashion  from  Scandinavia  to  Greenland  to  the  Pacific,  suggesting  that  they 
are  associated  with  Westward  propagating  planetary  waves.  When  the  study  was  repeated  using  a 
more  conventional  500  mb  zonal  wind  index,  the  results  did  not  show  such  a  clear  relationship  with 
the  teleconnection  patterns.  Thus,  Dr.  Weng  succeeded  in  developing  a  new  and  useful  index. 

2.  Baroclinic  Instability  in  the  Presence  of  Forcing 

The  focus  of  this  research  was  to  understand  the  interactions  between  the  stationary  planetary 
scale  and  the  transient  synoptic  scales  of  atmospheric  behavior.  To  this  end  Drs.  A.  Barcilon  and 
T.  Nathan  ( Geophya .  Astrophys.  Fluid  Dyn.,  56,  59-79,  1991)  studied  the  nonlinear  evolution  of 
a  free  baroclinic  wave  in  the  presence  of  slightly  supercritical,  vertically  sheared  zonal  flow  and  a 
forced  stationary  wave  field  consisting  of  a  single  zonal  scale  and  an  arbitrary  number  of  meridional 
harmonics.  The  analysis  was  carried  out  within  the  context  of  the  two-layer  baroclinic  model. 
Using  a  multiple  time  scale  analysis,  Drs.  Barcilon  and  Nathan  found  that  the  presence  of  the 
planetary-scale  stationary  wave  induces  zonal  variations  in  the  supercriticality  which  alters  the 
growth  rate  and  asymptotic  equlibration  of  the  baroclinic  wave.  This  comes  about  in  part  by  the 
direct  interaction  of  the  synoptic  scale  wave  with  the  forced  wave,  and  in  part  by  the  interaction 
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of  the  synoptic-scale  wave  with  the  zonal  mean  flow  which  has  been  altered  by  the  stationary  wave 
fluxes.  When  the  stationary  wave  field  is  confined  to  the  upper  layer  and  consists  of  only  the  gravest 
cross-stream  mode,  this  condition  is  favorable  for  the  non-zero  equilibration  of  the  free  wave.  When 
it  is  confined  to  the  lower  layer,  this  condition  is  unfavorable  for  such  equilibration. 

3.  Baroclinic  Wave  Structures  and  Transitions 

This  study  was  motivated  by  the  need  to  resolve  controversies  in  the  meteorological  literature 
concerning  the  effects  of  Ekman  layer  dissipation  and  sloping  boundaries  on  baroclinic  waves. 
For  example,  ( i )  Is  a  small  amount  of  Ekman  layer  dissipation  stabilize  or  destabilize  baroclinic 
waves?  (it)  If  there  is  frictional  instability  for  small  dissipation,  how  small  should  the  dissipation 
be?  (iii)  What  are  the  differences  between  the  effects  of  symmetric  and  asymmetric  Ekman  layer 
dissipation?  (tv)  What  are  the  effects  of  traditional  (more  friction  at  the  bottom  than  the  top) 
and  nontraditional  (less  friction  at  the  bottom  than  at  the  top  —  for  example,  if  there  is  strong 
gravity  waves  action  at  the  jet  stream  level)?  (v)  How  are  these  effects  altered  in  the  presence  or 
absence  of  sloping  boundaries?  (vi)  For  easterly  wind  shear  how  do  the  stability  results  change?  In 
an  effort  to  answer  these  questions,  Drs.  Weng  and  Barcilon  ( Ceophys .  Astrophys.  Fluid  Dyn .,  59, 
1-24,  1991)  first  studied  the  effects  of  Ekman  dissipation  and  sloping  boundaries  separately  in  a 
linear  framework  adding  factors  one-by-one  to  determine  how  each  altered  the  solution.  Then  they 
examined  the  effects  of  nonlinear  wave-mean  flow  and  wave-wave  interactions  on  baroclinic  waves 
in  the  presence  of  viscous  asymmetry  (preprint  of  Eighth  Conference  on  Atmospheric  and  Ocean 
Waves  and  Stability,  Oct.  14-18,  1991,  134-137).  Their  main  conclusions  were  as  follows: 

In  the  absence  of  sloping  boundaries,  symmetric  Ekman  layer  dissipation  has  a  stabilizing 
tendency  for  all  waves.  Asymmetric  Ekman  layer  dissipation,  on  the  other  hand,  destabilizes  the 
waves  by  extending  the  unstable  wave  band  toward  both  long  and  short  waves  and  by  increasing 
all  growth  rates  in  comparison  with  symmetric  dissipation.  It  also  makes  the  waves  dispersive  and 
their  structure  asymmetric  about  mid-depth. 

In  the  presence  of  sloping  boundaries  wave  dispersion  increased  and  the  wave  propagation 
characteristics  became  sensitive  to  the  sign  of  the  shear  and  the  sense  of  the  viscous  asymmetry. 
For  westerly  wind  shear,  the  wave  amplitude  was  reduced  and  the  waves  became  more  baroclinic  in 
the  lower  atmosphere  and  barotropic  above.  The  reverse  wave  true  for  easterly  wind  shear.  The 
inclusion  of  nonlinearity  did  not  affect  these  basic  results. 
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4.  Role  of  Mountains 


The  goal  of  this  research  was  to  ascertain  the  effects  of  topographically  forced  planetary  waves 
on  the  genesis  and  evolution  of  mid-latitude  synoptic  disturbances.  The  approach  was  to  inves¬ 
tigate  the  linear  and  nonlinear  stability  of  zonally  varying  basic  flows  within  the  context  of  the 
quasi-geostrophic  barotropic  and  two-layer  baroclinic  models.  The  models  were  constructed  so  as 
to  combine  baroclinic,  barotropic  and  topographic  processes  in  the  form  of  primary  and  secondary 
instabilities.  The  treatment  of  the  nonlinear  problem  was  done  analytically  using  asymptotic  ex¬ 
pansions  and  multiple  time  scales. 

To  this  end,  Drs.  Barcilon  and  Nathan  first  examined  the  nonlinear  evolution  of  a  free  baro¬ 
clinic  wave  in  a  two-layer  atmosphere  in  the  presence  of  slightly  supercritical,  vertically  sheared 
zonal  flow  and  a  forced  stationary  wave  consisting  of  a  single  zonal  scale  with  an  arbitrary  number 
of  meridional  harmonics.  Using  a  multiple  time  scale  analysis,  they  showed  that  the  stationary 
wave  induces  zonal  variations  in  supercriticality  which  alters  the  growth  rate  and  asymptotic  equi¬ 
libration  of  synoptic  disturbances  via  two  mechanisms.  The  first  is  the  direct  interaction  of  the 
synoptic  disturbance  with  the  forced  wave  (i.e.,  a  wave-wave  mechanism).  The  second  is  the  in¬ 
teraction  of  the  synoptic  disturbance  with  the  portion  of  the  mean  field  that  has  been  changed  by 
the  stationary  wave  fluxes  (a  wave-mean  flow  mechanism).  These  mechanisms  oppose  or  augment 
each  other  depending  on  the  amplitude  and  spatial  structure  of  the  stationary  wave  field.  If  it 
is  confined  to  the  lower  layer,  and  consists  only  of  the  gravest  cross  stream  mode,  the  free  wave 
does  not  equilibrate  to  a  finite  amplitude.  If  it  is  confined  to  the  upper  layer,  the  free  wave  does 
equilibrate  to  a  finite  amplitude.  The  interaction  of  a  synoptic  wave  with  a  stationary  wave  con¬ 
sisting  of  two  or  more  meridional  harmonics  generates  time  dependent  heat  fluxes  that  vary  with 
the  period  of  the  synoptic  wave.  If,  however,  the  stationary  wave  field  contains  several  sufficiently 
large  amplitude  meridional  harmonics  it  can  destroy  the  synoptic  wave. 

In  a  further  study  ( Trend s  in  Atmos.  Sci .,  1992,  in  press),  Dr.  Nathan  addressed  the  issue  of 
how  the  combined  influence  of  topographic  and  zonally  varying  potential  vorticity  (PV)  forcing  can 
influence  low  frequency  motions  in  the  atmosphere.  In  particular,  he  demonstrated  analytically  that 
such  forcing  can  produce  dramatic  changes  in  the  baroclinic  topographic  instability  properties  of 
the  flow.  The  PV  forcing  can  enhance,  suppress  or  catalyze  the  topographic  instability  depending 
on  the  phase  of  the  PV  fordng  with  respect  to  the  topography.  These  changes  result  from  the 
alteration  of  the  zonal  mean  flow  produced  by  the  interaction  between  the  resonant  wave  and 
the  PV  fordng.  Since  the  PV  fordng  in  Dr.  Nathan’s  model  reflects  the  diabatic  heating  in  the 
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atmosphere,  a  major  conclusion  of  this  study  is  that  land-  sea  heating  contrasts  may  play  a  more 
important  role  in  the  dynamics  of  topographically  induced  instabilities  than  previously  thought. 

Drs.  Nathan  and  Barcilon  also  extended  the  Jin  and  Ghil  (1990)  barotropic  model  of  topo¬ 
graphically  induced  instabilities  to  include  a  steady,  zonaliy  varying  vorticity  source.  In  particular, 
they  showed  that,  in  contrast  to  the  results  obtained  by  JG,  the  vorticity  forcing  is  capable  of 
generating  low  frequency  oscillations  that  remain  strictly  upstream  of  the  topographic  ridge.  For 
sufficiently  large  vorticity  forcing,  the  Hopf  bifurcation  is  destroyed,  irrespective  of  the  relative 
phase  of  the  forcings.  In  addition  to  the  lower  frequency  oscillations  obtained  by  JG,  they  also 
obtained  solutions  that  exhibit  two  distinctly  different  types  of  regimes.  One  regime  is  a  westward- 
propagating,  low  frequency  oscillation  that  has  the  character  of  a  free  wave.  The  second  regime  is 
characterized  by  a  quasi-steady  state  in  which  the  wave  ridge  remains  upstream  of  the  topographic 
ridge.  This  two  regime  solution  is  quite  robust  and  only  occurs  when  the  vorticity  forcing  augments 
the  topographic  forcing. 

Finally,  Drs.  Barcilon  and  Blumen  ( Geophys .  Astrophys.  Fluid  Dyn.,  58,  25-43, 1991)  exam¬ 
ined  the  formation  of  large-amplitude  downslope  wind  storms  such  as  those  that  occur  to  the  lee 
of  the  Rockies.  They  found  criteria  on  the  vertical  stratification  of  the  atmosphere  and  cross-ridge 
wind  variations  that  would  cause  internal  gravity  waves  generated  by  an  airstream  passing  over  a 
mountain  to  amplify  explosively.  The  cross-ridge  wind  variations  produce  significant  differences  in 
the  surface  wind  in  the  lee  of  the  mountain  and  in  the  surface  drag  compared  to  cases  with  no  such 
variations. 

5.  Complex  Principal  Component  Analysis 

Low  frequency  oscillations  were  also  studied  by  Dr.  Pfeffer  and  his  colleagues  (Bernardet,  et. 
al.,  J.  Atmos.  Sci..,  47,  3023-3043,  1990)  using  complex  principal  component  analysis  of  data 
from  rotating  laboratory  experiments  in  a  thermally  driven  fluid  over  bottom  topography.  In  these 
experiments  the  topography  generated  a  planetary-scale  stationary  wave  and  the  barodinicity  pro¬ 
duced  synoptic-scale  traveling  disturbances.  The  complex  principal  component  analysis  of  data 
measured  at  two  levels  in  the  fluid  showed  the  presence  of  the  standing  wave  with  a  baroclinic 
structure  which  modulated  the  synoptic  scale  waves  traveling  over  it.  The  analysis  also  revealed 
that  the  travelling  wave  energy  was  concentrated  downstream  of  the  topographic  ridges,  suggesting 
the  mechanism  by  which  storm  tracks  are  maintained  in  the  atmosphere  —  primarily  by  mountains 
rather  than  by  land-sea  temperature  contrasts.  Another  significant  result  was  that  of  symmetry 
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breaking.  That  is,  although  the  topography  consisted  of  two  perfectly  symmetric,  sinusoidal  moun¬ 
tains  the  wave  trains  downstream  of  one  mountain  were  quite  different  from  those  downstream  of 
the  other  suggesting  that  it  is  not.  The  different  shapes  of  the  Rockies  and  the  Himalayas  that 
decouple  persistent  anomalies  in  the  Atlantic  from  those  in  the  Pacific,  but  rather  the  nonlinear 
dynamics  that  lead  to  symmetry  breaking  even  in  simple  laboratory  experiments. 

The  data  from  another  series  of  laboratory  experiments  was  also  used  by  Drs.  Eisner  and 
Tsonis  to  test  methodologies  of  nonlinear  prediction  of  chaotic  time  series. 

6.  Numerical  Studies 

Our  numerical  studies  were  motivated  by  the  need  to  supplement  nonlinear  analytic  stud¬ 
ies  with  numerical  studies  and  to  understand  the  limitations  of  both.  In  this  connection  Dr. 
Navon  and  his  student  Jian  Yu  ( Computers  in  Geosciences ,  17,  1311-1343,  1991)  developed  and 
documented  a  computer  program  applying  the  Turkel-Zwas  explicit  large  time-step  scheme  to  a 
hemispheric  barotropic  model  with  constraint  restoration  of  integral  invariants  of  the  shallow  water 
equations.  Mr.  Yu  (FSU-SCRI-90T-124,  1990)  also  implemented  and  documented  the  software 
package  PLTMG  for  solving  a  general  second  order  parameter  involved  system  of  elliptic  partial 
differential  equations  on  a  two-dimensional  domain  with  Dirichlet  and/or  Neuman  boundary  condi¬ 
tions.  And  Drs.  Pfeffer,  Navon  and  Zou  studied  the  effects  of  using  two  different  time-differencing 
schemes  in  the  numerical  integration  of  a  global  model.  This  last  study  revealed  that  significant 
differences  in  the  solutions  occur  when  only  the  time  scheme  is  changed.  This  makes  it  imperative 
to  carry  out  both  analytical  and  numerical  studies  and  to  limit  conclusions  only  to  results  that  can 
be  corroborated  by  both  methodologies. 
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III.  ABSTRACTS  OF  PAPERS 
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Contribution  No  "2-  ^ 
Gacphysicd  Ftuid  Dy'^TbTtitute 
i  lorioa  biate  University 

Reprinted  from  Journal  of  the  atmospheric  Sciences,  Vol.  47.  No.  24, 1 5  December  1 990 

Amoiciti  McteorotogfOi  Socxi* 


Low-Frequency  Oscillations  in  a  Rotating  Annulus  with  Topography 

P.  Bernardet,  a.  Butet  and  M.  D£qu£ 

Centre  National  de  Recherches  Meteoroiogiques,  EERM/DMN.  Toulouse.  France 


M.  GHIL 

Climate  Dynamics  Center.  Department  of  Atmospheric  Sciences  and  Institute  of  Geophysics  and  Planetary  Physics. 

University  of  California  at  Los  Angeles.  California 


R.  L.  PFEFFER 

Geophysical  Fluid  Dynamics  Institute  and  Department  of  Meteorology.  Florida  State  University.  Tallahassee.  Florida 
(Manuscript  received  19  February  1990,  in  final  form  13  July  1990) 

ABSTRACT 

Experiments  were  performed  in  a  rotating,  differentially  heated  annulus,  with  and  without  bottom  topography 
of  azimuthal  wavenumber  2.  Both  water  and  a  viscous  glvcerol-water  mixture  were  used  as  a  working  fluid,  (n 
one  series  of  experiments,  measurements  of  azimuthal  velocity  u  were  carried  out  by  Doppler-laser  velocimetry 
at  midradius  and  at '/j  and  :/j  depth.  In  the  other,  temperature  measurements  were  made  by  a  set  of  thermistors 
at  three  different  heights  and  three  different  radii.  Results  were  analyzed  by  Fourier  transformation,  separatelv 
in  space  and  in  time,  and  in  terms  of  complex  empirical  orthogonal  functions  (CEOFs). 

In  the  experiments  with  topography,  a  standing  wave  2  is  generated,  with  larger  amplitude  at  the  upper  level 
and  a  tilted  wave  structure.  The  two  leading  CEOFs  contain  a  very  large  fraction  of  the  variance,  and  give  an 
excellent  picture  of  the  spatial  modulation  of  the  traveling  baroclinic  waves.  The  dominant  barodinic  wave  has 
azimuthal  wavenumber  4.  5  or  6.  according  to  the  nondimensional  parameters  of  the  given  experiment,  and 
pronounced  sidebands  due  to  the  topography.  The  modulation  of  this  wave  is  such  that  its  largest  amplitude 
occurs  at  the  lower  level  upstream  of  the  two  topographic  ridges.  At  the  upper  level,  the  modulation  is  weaker, 
with  the  maximum  wave  amplitude  located  downstream  of  the  ridges.  Partial  decoupling  of  the  two  wave  trains 
attached  to  the  two  ridges  is  evident  in  one  experiment 

Low-frequencv  vacillation  of  the  entire  flow  pattern  is  apparent:  this  vacillation  has  a  period  of  about  $0 
annulus  rotations  in  the  viscous  mixture.  The  possible  relevance  of  this  topographically  induced  vacillation  to 
the  extratropical  30-60  day  oscillation  is  discussed. 


1.  Introduction  and  motivation 

Meteorology  is  mainly  an  observational,  rather  than 
an  experimental  science.  Still,  rotating  annulus  exper¬ 
iments  represent  a  powerful  paradigm  for  the  general 
circulation  of  the  atmosphere’s  midlatitudes,  including 
its  variability  (Lorenz  1967,  Ch.  VI;  Ghil  and  Childress 
1 987,  Ch.  5 ) .  Detailed  comparisons  between  the  fluid- 
flow  phenomena  and  their  causal  mechanisms  in  the 
annulus,  on  the  one  band,  and  analogous  phenomena 
in  the  atmosphere,  on  the  other,  are  often  difficult  and 
sometimes  controversial  (Wallace  and  Hsu  1985).  But 
the  rotating,  differentially  heated  annulus  provides  a 
fluid  environment  for  which  parameters  can  be 
changed,  measurements  taken  and  observed  phenom¬ 
ena  explained  with  greater  ease  than  for  the  atmosphere 
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itself.  Whether  these  phenomena,  and  their  explana¬ 
tions,  are  the  same  as  in  the  atmosphere  has  to  be  de¬ 
cided  on  a  case-by-case  basis.  The  dynamically  inter¬ 
esting  observations  and  explanations,  however,  arising 
from  the  annulus  experiments  are  an  important  source 
of  ideas  in  the  study  of  large-scale  atmospheric  motions. 

In  studies  of  the  atmosphere,  low-frequency  vari¬ 
ability  occupies  an  increasingly  important  part  (Hos¬ 
kins  and  Pearce  1983;  Branstator  et  al.  1988).  In  par¬ 
ticular,  intraseasonal  oscillations  with  periods  of 
roughly  25  days  (Branstator  1987;  Ghil  and  Mo  1990) 
and  40  days  ( Weickmann  et  aL  1985;  Lau  and  Phillips 
1986;  Ghil  and  Mo  1990)  in  the  Northern  Hemisphere 
(NH)  extratropics  have  attracted  some  attention  re¬ 
cently.  Theoretical  explanations  of  the  extratropical 
NH  30-60  day  oscillation  have  tended  to  fall  into  two 
categories.  First,  these  oscillations  have  been  seen  as 
an  essentially  linear  response  to  tropical  forcing  with 
the  same  range  of  periods  (see,  for  example.  Hsu  et  al. 
1990,  and  references  therein).  Second,  they  have  been 
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The  weakly  nonlinear  evolution  of  a  free  baroclinic  wave  in  the  presence  of  slightly  supercritical, 
vertically  sheared  zonal  flow  and  a  forced  stationary  wave  field  that  consists  of  a  single  zonal  scale  and 
an  arbitrary  number  of  meridional  harmonics  is  examined  within  the  context  of  the  conventional  two- 
layer  model.  The  presence  of  the  (planetary-scale)  stationary  wave  introduces  zonal  variations  in  the 
supercriticality  and  is  shown  to  alter  the  growth  rate  and  asymptotic  equilibrium  of  the  (synoptic-scale) 
baroclinic  wave  via  two  distinct  mechanisms:  The  first  is  due  to  the  direct  interaction  of  the  stationary 
wave  with  the  shorter  synoptic  wave  (wave-wave  mechanism),  and  the  second  is  due  to  the  interaction 
of  the  synoptic  wave  with  that  portion  of  the  mean  Reid  that  is  corrected  by  the  zonally  rectified 
stationary  wave  fluxes  (wave-mean  mechanism).  These  mechanisms  can  oppose  or  augment  each  other 
depending  on  the  amplitude  and  spatial  structure  of  the  stationary  wave  field.  If  the  stationary  wave 
field  is  confined  primarily  to  the  upper  (lower)  layer  and  consists  of  only  the  gravest  cross-stream 
mode,  conditions  are  favorable  (unfavorable)  for  nonzero  equilibrium  of  the  Tree  wave. 

In  addition  to  the  time  dependent  heat  flux  generated  by  baroclinic  growth  of  the  free  wave,  its 
interaction  with  a  stationary  wave  field  consisting  of  two  or  more  meridional  harmonics  generates  time 
dependent  heal  fluxes  that  vary  with  period  of  the  Tree  wave.  However,  if  the  stationary  wave  field 
contains  several  meridional  harmonics  of  sufficiently  large  amplitude,  the  free  baroclinic  wave  is 
destroyed. 

KEY  WORDS:  Barodinic  instability,  stationary  waves,  wave-wave  and  wave-mean  flow  interactions. 
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1.  INTRODUCTION 

A  characteristic  signature  of  geophysical  fluid  systems  is  the  presence  of  several 
mutually  interacting  scales  of  motion  (e.g^  Colucci  tt  aU  1981;  Li  1986; 

Pfeffer  et  al ,  1990).  Here  we  focus  on  the  finite-amplitude  interactions  between  two 
scales  of  motion  that  characterize  the  large-scale  circulation  of  the  atmosphere:  the 
stationary  planetary-scales  and  the  transient  synoptic-scales. 
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ASYMMETRIC  EKMAN  DISSIPATION, 
SLOPING  BOUNDARIES  AND 
LINEAR  BAROCLINIC  INSTABILITY 
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The  effects  of  symmetric  and  asymmetric  Ekman  dissipation  on  baroclinic  instability,  phase  speed  and 
wave  structure  in  a  linear  Eady-like  channel  model  with/without  oppositely  sloping  boundaries  are 
examined.  In  the  absence  of  sloping  boundaries,  symmetric  Ekman  dissipation  has  a  stabilizing  tendency 
for  all  waves.  When  the  Ekman  layers  are  asymmetric,  the  viscous  asymmetry  destabilizes  waves  by  (i) 
extending  the  unstable  waveband  toward  both  long  and  short  waves,  and  (ii)  increasing  their  growth 
rates,  compared  with  viscous  symmetric  case.  When  the  asymmetric  dissipation  is  very  small,  the 
destabilization  may  result  in  a  frictional  instability  for  short  waves  which  are  stable  in  inviscid  case.  The 
viscous  asymmetry  makes  waves  dispersive  and  their  structure  asymmetric  about  mid-depth.  However, 
the  sense  of  the  viscous  asymmetry  and  the  sign  of  shear  do  not  affect  the  instability,  but  do  modify  the 
direction  of  phase  propagation  and  the  shape  of  the  wave  structure. 

With  asymmetric  Ekman  dissipation  and  sloping  boundaries,  frictional  instability  is  a  mixture  of  three 
mechanisms:  (i)  symmetric  dissipation  in  the  presence  of  sloping  boundaries:  (ii)  viscous  asymmetry  in 
the  absence  of  sloping  boundaries:  and  (iii)  the  sense  of  viscous  asymmetry  in  the  presence  of  sloping 
boundaries,  which  is  sensitive  to  the  wavenumber  and  the  sign  of  shear.  The  slope  renders  the  phase  speed 
more  dispersive  and  sensitive  to  the  sign  of  shear  and  the  sense  of  the  viscous  asymme^v.  For  westerly 
shear,  the  slope  reduces  (increases)  the  wave  amplitude  and  makes  the  wave  more  baroclinic  (barotropic) 
in  the  lower  (upper)  level.  For  a  given  wavenumber,  easterly  shear  dynamics  may  be  deduced  from  westerly 
shear  dynamics  by  a  proper  change  of  the  sense  of  the  viscous  asymmetry. 

KEYWORDS:  Baroclinic  waves,  Ekman  dissipation,  sloping  boundaries. 

1.  INTRODUCTION 

Because  of  its  linear  and  nonlinear  analytical  simplicity,  the  Eady  (1949)  model  of 
baroclinic  instability  has  served  as  a  "workhorse"  to  test  various  aspects  of  that 
instability.  Yet,  we  feel  that  the  various  possible  alternatives  offered  by  this  model, 
such  as  asymmetric  Ekman  dissipation  and  sloping  boundaries  in  particular,  could 
stand  a  thorough  re-examination  which  we  offer  below  to  lay  the  foundation  before 
proceeding  with  nonlinear  studies. 

Using  a  two-layer  /f-plane  model,  Holopainen  (1961)  considered  a  single  Ekman 
layer  at  the  lower  boundary  and  found  that  the  effect  of  friction  reduces  the  amplitude 
of  the  perturbation,  except  within  two  narrow  wavebands  adjacent  to  the  two  inviscid 
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The  steady,  hydrostatic,  inviscid,  Boussinesq  flow  of  a  stably  stratified  fluid  over  a  bell-shaped  ridge  is 
investigated  within  the  framework  of  a  linear  model  The  three  layer  model  atmosphere  introduced  is 
such  that  the  Brunt- Viisali  frequency  is  constant  in  each  layer  but  the  interfaces  of  the  middle  layer 
are  allowed  to  vary  gently  in  the  cross-ridge  direction.  In  essence,  the  model  can  be  tuned  in  both 
vertical  and  horizontal  directions.  These  cross-ridge  variations  can  produce  significant  differences  in 
both  the  cross-ridge  surface  wind  and  the  surface  drag  compared  to  the  response  obtained  by  use  of  a 
horizontally  uniform  reflecting  layer.  These  changes  are  sensitive  to  both  the  vertical  location  of  the 
middle  layer  and  to  the  slope  of  its  lower  interface  at  the  ridge  crest.  Many  or  these  features  are 
explained  by  means  of  a  conventional  layered-model  analysis. 

KEY  WORDS:  Mountain  waves,  stratified  flow,  Boussinesq  fluid. 

1.  INTRODUCTION 

Typically,  the  static  stability,  expressed  by  the  Brunt-Vaisala  frequency,  undergoes 
considerable  variability  in  the  vertical  direction.  Figure  I  shows  that  the  static 
stability,  N2,  exhibits  variability  over  scales  of  1-2  kin  in  the  troposphere.  Similar 
features  are  retained  between  the  two  stations,  as  in  the  region  between  6  and 
8  km  in  panel  b.  In  other  cases,  however,  the  apparent  layering  of  the  static 
stability  may  be  localized  in  the  vicinity  of  each  station,  and/or  the  same  features 
may  ascend  or  descend  in  altitude  between  GCT  and  DEN.  This  latter  possibility, 
which  may  characterize  the  features  in  panel  a,  is  examined  in  the  present  study. 

We  consider  the  dynamics  of  hydrostatic  mountain  waves  in  the  presence  of 
static  stability  variations  in  the  cross-ridge  direction.  Such  states  will  consist  of 
regions  bounded  by  x-dependent  interfaces  across  which  the  Scorer  parameter 
changes  abruptly,  remaining  constant  on  either  sides  of  these  interfaces.  These 
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A  potential  vorticity  (PV)  index  is  defined  as  a  measure  of  the  zonally  averaged,  mid-latitude  PV  gradient 
on  the  300  K  isentropic  surface  in  the  northern  hemisphere,  and  use  is  made  of  the  1978/79  winter  ECMWF 
FGGE  m-b  data  set  to  study  the  evolution  of  that  index  and  its  relation  to  teleconnection  patterns  of  500  mb 
geopotential-height  anomaly.  This  index  vacillates  with  a  dominant  period  of  about  two  weeks.  For  some  time 
lags,  the  isoline  patterns  of  the  cross-correlation  function  between  the  PV  index  and  500  mb  geopotential- 
height  anomalies  are  similar  to  some  teleconnection  patterns  studied  by  several  investigators  using  large 
climatological  data  sets.  Scandinavia,  Greenland  and  the  Pacific  are  important  centres  of  ‘action’  in  the  northern 
hemisphere  in  the  1978/79  winter.  The  formation  of  blocking  highs  at  these  locations  occurs  at  rather  specific 
phases  during  the  vacillatory  evolution  of  the  PV  index,  especially  as  the  flow  transits  from  a  low  to  a  high  PV 
index. 

When  the  above  study  is  duplicated  using  a  zonal  mean  wind  gradient  at  500  mb  as  a  zonal  index,  the 
results  do  not  reflect  the  teleconnection  patterns  as  crisply  as  the  PV  index  does  for  the  FGGE  year. 


1.  Introduction 

Using  the  European  Centre  for  Medium-range  Weather  Forecast  (ECMWF)  First 
GARP  Global  Experiment  (FGGE)  III-b  data  from  1  February  to  15  May  1979,  Kidson 
(1985)  defined  a  zonal  index  in  terms  of  the  zonal  mean  wind  difference  between  two 
latitude  belts  at  500  mb,  one  centred  at  60®N  and  the  other  centred  at  39.25“N.  His 
spectral  analysis  showed  that  the  20-  to  30-day  cycle  was  most  prominent  near  40°N- 
60°N.  He  found  that  a  low-index  flow  related  to  the  presence  of  persistent  cut-off  lows 
at  500  mb  off  the  coast  of  Newfoundland,  with  a  warm  ridge  or  a  cut-off  high  in  the 
vicinity  of  Iceland;  while  a  high-index  flow  corresponded  to  a  zonal  flow  at  high  latitudes, 
with  slow-moving  patterns  over  the  Atlantic  at  lower  latitudes  and  to  the  south  of  the 
main  jet.  Although  his  correlation  pattern  between  the  mean-sea-level  pressure  and  the 
zonal  index  was  broadly  similar  to  the  North  Atlantic  Oscillation  (NAO)  pattern  found 
in  climate  data,  some  differences  were  apparent,  particularly  the  absence  of  any  significant 
variations  in  the  Pacific,  which,  if  present,  should  be  related  to  the  North  Pacific 
Oscillation  (NPO)  or/and  the  Pacific/North  American  (PNA)  patterns. 

Teleconnections  and  seesaws  have  been  found  in  climate  studies  by  many  inves¬ 
tigators.  For  example,  van  Loon  and  Rogers  (1978,  hereafter  referred  to  as  vLR)  found 
the  Greenland-Northern  Europe  seesaw  pattern,  where  winter  temperatures  are  below 
normal  over  northern  Europe  while  being  above  normal  over  Greenland  and  the 
Canadian  arctic;  and  pice  versa.  They  also  found  that  sea-level-pressure  anomalies  over 
most  of  the  northern  hemisphere  are  associated  with  this  regional  seesaw  in  temperature. 
Wallace  and  Gutzler  (1981,  hereafter  referred  to  as  WG)  summarized,  based  on  existing 
literature,  four  teleconnection  patterns;  the  NAO,  the  NPO,  the  PNA.  and  a  zonally 
symmetric  seesaw  pattern  found  by  Lorenz  (1951)  in  the  zonally  averaged  sea-level- 
pressure  field.  Since  the  strength  and  location  of  these  patterns  vary  from  winter  to 
winter,  the  main  climatic  teleconnection  patterns  may  differ,  depending  on  the  data  set 
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EXSHALL:  A  TURKEL-ZWAS  EXPLICIT  LARGE 
TIME-STEP  FORTRAN  PROGRAM  FOR  SOLVING 
THE  SHALLOW- WATER  EQUATIONS 
IN  SPHERICAL  COORDINATES* 
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Abstract — A  FORTRAN  computer  program  is  presented  and  documented  applying  the  Turkel-Zwas 
explicit  large  time-step  scheme  to  a  hemispheric  barotropic  model  with  constraint  restoration  of  integral 
invariants  of  the  shallow-water  equations.  We  then  proceed  to  detail  the  algorithms  embodied  in  the  code 
EXSHALL  in  this  paper,  particularly  algorithms  related  to  the  efficiency  and  stability  of  T-Z  scheme  and 
the  quadratic  constraint  restoration  method  which  is  based  on  a  variational  approach.  In  particular  we 
provide  details  about  the  high-latitude  filtering,  Shapiro  filtering,  and  Robert  filtering  algorithms  used  in 
the  code.  We  explain  in  detail  the  various  subroutines  in  the  EXSHALL  code  with  emphasis  on  algorithms 
implemented  in  the  code  and  present  the  flowcharts  of  some  major  subroutines.  Finally,  we  provide  a 
visual  example  illustrating  a  4-day  run  using  real  initial  data,  along  with  a  sample  printout  and  graphic 
isoline  contours  of  the  height  field  and  velocity  fields. 

Key  Words :  Shallow-water  equations.  Spherical  coordinates.  Explicit  finite-differences.  Constraint  restor¬ 
ation,  Filtering  techniques. 


INTRODUCTION 

Recently,  a  considerable  amount  of  work  has  -been 
dedicated  and  aimed  at  efficient  integration  of  shal¬ 
low-water  equations  in  view  of  using  these  methods 
in  numerical  weather  prediction  models.  In  order 
to  achieve  computational  accuracy  and  efficiency, 
most  methods  are  concerned  with  the  different  time- 
scale  of  the  advection  and  the  gravity-inertia  terms 
in  the  shallow-water  equations  model  separately. 
Semi-implicit  schemes  (Robert,  1979;  Burridge,  1975) 
and  split-explicit  schemes  (Magazenkov,  Shvets, 
and  Shneyerov,  1971;  Gadd,  1978a,  1978b)  are 
examples  of  those  methods.  In  the  split-explicit 
schemes,  a  substantial  computational  economy  is 
achieved  when  compared  to  usual  explicit  time 
integration  schemes. 

Turkd  and  Zwas  (1979)  proposed  a  space-splitting 
rather  than  a  time-splitting  algorithm  for  the  explicit 
integration  of  the  shallow-water  equations.  Their 
method  is  based  on  the  fact  that  the  fast  gravity- 
inertia  waves  contain  only  a  small  fraction  of  the 
total  available  energy  and  therefore  these  waves  can 
be  calculated  with  a  lower  accuracy  than  the  slow 
Rossby  waves,  that  is  on  a  coarser  mesh.  An  appli¬ 
cation  of  the  T-Z  space  split-explicit  integration 
schemes  with  real  initial  data  is  presented  and  dis- 
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cussed  by  Navon  and  de  Villiers  (i987)  and  its 
properties  are  discussed  further  in  Neta  and  Navon 
(1989).  A  linear  transfer  function  analysis  of  the 
shallow-water  equations  in  spherical  coordinates  for 
the  Turkel-Zwas  explicit  large  time-step  scheme  was 
carried  out  by  Neta.  Navon.  and  Yu  (1990). 

The  purpose  of  this  paper  is  to  present  a  practical 
FORTRAN  code,  EXSHALL,  which  implements 
the  T-Z  explicit  large  time-step  scheme  for  the  shal¬ 
low-water  equations  in  spherical  coordinates  along 
with  constraint  restoration  methods  for  enforcing  a 
posteriori  conservation  of  the  integral  invariants  of 
the  shallow-water  equations.  The  computer  program 
is  explained  in  detail  in  connection  with  the  various 
algorithms  implemented  in  the  code  EXSHALL. 
This  paper  can  be  used  as  a  user’s  guide  to  the 
program  EXSHALL  both  in  providing  a  brief 
description  of  the  theory  as  well  as  detailed  program¬ 
ming  implementation. 

We  present  here  the  T-Z  scheme  for  the  shallow- 
water  equations  in  spherical  coordinates  and  its 
related  algorithmic  background.  Various  filters  used 
with  T-Z  scheme  and  which  impact  on  its  stability 
also  are  presented;  a  detailed  description  of  the 
various  subroutines  in  the  code  EXSHALL  is  pre¬ 
sented;  and  a  typical  example  of  a  4-day  run  with  the 
program  EXSHALL  is  presented  along  with  graphi¬ 
cal  output  Finally,  in  the  Appendix  the  commented 
and  documented  FORTRAN  source  listing  the  code 
of  the  program  EXSHALL  is  attached. 
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ABSTRACT 

This  paper  presents  evidence  of  the  sensitivity  of  a  general  circulation  model  (GCM)to  the  time-differencing 
scheme  employed  when  the  physical  parameterization*  and  space  discretization  are  not  changed.  For  this  purpose, 
two  time-marching  schemes — the  leapfrog  and  the  Matsuno  schemes — are  analyzed  and  tested  on  the  National 
Aeronautics  and  Space  Administration-Goddard  Laboratory  for  Atmospheric  Studies  (NASA-GLAS)  fourth- 
order  GCM  in  terms  of  the  stability  and  behavior  of  2-month-averaged  fields.  Linear  analysis  suggests  that 
Rossby  waves  are  slightly  damped  and  slightly  accelerated  when  the  Matsuno  scheme  is  used  and  that  these 
effects  are  scale  selective,  being  smallest  for  the  longest  waves.  It  also  suggests  that  such  waves  are  accelerated 
less  and  are  not  damped  when  the  leapfrog  scheme  is  used.  An  empirical  orthogonal  function  analysis  of  the 
meridional  component  of  velocity  at  46*N,  keeping  at  least  70%  of  the  variance,  reveals  less  shortwave  activity 
in  the  numerical  solution  with  the  Matsuno  scheme  but  does  not  lend  support  to  the  conclusion  that  the  waves 
are  accelerated  less  in  the  solution  with  the  leapfrog  scheme. 

The  two-dimensional  EBasscn-Palm  (E-P)  flux  divergence  and  the  eddy-induced  mean  meridional  diculanon 
are  found  to  be  stronger  in  the  simulation  with  the  leapfrog  rime-differencing  scheme  than  in  the  one  with  the 
Matsuno  scheme,  suggesting  that  the  transient-wave  activity  is  damped  by  the  Matsuno  scheme.  On  the  other 
hand,  the  three-dimensional  starionary-wave  activity  flux  in  the  Northern  Hemisphere  simulated  with  the  Matsuno 
scheme  is  more  intense  than  that  simulated  with  the  leapfrog  scheme,  indicating  that  the  stationary  waves  are 
more  robust  in  the  integration  with  the  Matsuno  scheme. 

The  GCM  precipitation  when  integrated  with  the  leapfrog  scheme  is  much  more  intense  over  the  tropical 
western  Pacific  and  the  northeastern  Pacific  and  less  intense  over  the  western  North  Atlantic  Ocean.  Tbc  kinetic 
energy  of  waves  with  wavenumber  greater  than  9  simulated  by  the  Matsuno  scheme  is  consistently  smaller  than 
that  obtained  by  the  leapfrog  scheme.  These  results  give  evidence  that  climate  simulations  are  sensitive  not  only 
to  physical  parameterizatioas  of  subgrid-scale  processes  but  also  to  the  numeral  methodology  employed. 


1.  Introduction 

General  circulation  models  (GCMs)  are  the  most 
elaborate  of  a  hierarchy  of  mathematical  models  used 
in  the  study  of  climate.  A  genera!  review  of  atmospheric 
GCMs  presented  by  Simmons  and  Bengtsson  ( 1984) 
emphasizes  the  importance  of  physical  processes  in  de¬ 
termining  the  behavior  of  climatic  systems.  GCMs  have 
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been  used  for  seasonal  simulation  (Shukla  et  al.  1981) 
as  well  as  for  studying  climate  variability  on  time  scales 
of  a  month  and  upward.  In  particular,  studies  have 
been  conducted  to  determine  the  sensitivity  of  atmo¬ 
spheric  GCMs  to  changes  in  physical  mechanisms  such 
as  surface  albedo  (Chamey  197S),  sea-ice  limits  in  the 
Arctic  (Herman  and  Johnson  1978),  low-frequency 
variability  (Chamey  and  Shukla  1981),  and  inadequate 
orographic  effects  (Wallace  et  aL  1983). 

Reviews  of  the  numerical  techniques  used  in  nu¬ 
merical  weather  prediction  models  and  GCMs  have 
been  presented  by  Mesinger  and  Arakawa  ( 1976)  and 
Kasahara  ( 1979).  Numerical  experiments  in  which  the 
earth’s  climate  is  simulated  using  different  numerical 
schemes  can  play  a  valuable  role  in  clarifying  the  nature 
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1.  INTRODUCTION 

The  effect  of  asymmetric  Ekman  dissipation  on  baro- 
dinic  waves  has  been  studied  in  linear  (e.g.,  Holopainen, 
1961)  and  nonlinear  (e.g.,  Mak,  1987)  frameworks.  Weng 
and  Barcilon  (1991,  VVB)  compared  the  effects  of  asymmet¬ 
ric  Ekman  dissipation,  with/without  sloping  boundaries,  on 
//near  barodinic  waves  in  an  Eady-like  model.  They  found 
that  when  the  asymmetric  Ekman  layers  are  horizontal,  vis¬ 
cous  asymmetry  destabilizes  waves  by  (i)  extending  the  un¬ 
stable  waveband  toward  both  long  and  short  waves,  and  (ii) 
increasing  their  growth  rates,  when  compared  with  the  sym¬ 
metric  case.  When  the  asymmetric  dissipation  is  very  small, 
the  destabilization  may  result  in  a  frictional  instability  for 
short  waves  which  are  stable  in  the  invisdd  case.  Viscous 
asymmetry  makes  waves  dispersive  and  their  structures  are 
asymmetric  about  mid-depth.  However,  the  sense  of  viscous 
asymmetry  and  the  sign  ot  the  vertical  shear  do  not  affect  the 
instability,  but  do  modify  the  direction  of  phase  propagation 
and  the  shape  of  the  wave  structure.  When  the  asymmet¬ 
ric  Ekman  layers  are  oppositely  sloping,  the  effect  of  viscous 
asymmetry  is  modified  by  the  slope.  The  frictional  instability 
is  a  mixture  of  three  mechanisms:  (i)  symmetric  dissipation 
in  the  presence  of  sloping  boundaries;  (ii)  viscous  asymme¬ 
try  in  the  absence  of  sloping  boundaries;  and  (iii)  the  sense 
of  viscous  asymmetry  in  the  presence  of  sloping  boundaries, 
which  is  sensitive  to  the  wavenumber  and  the  sign  of  the 
shear.  Thus,  the  slope  renders  the  phase  speed  more  disper¬ 
sive  and  sensitive  to  the  sign  of  the  shear  and  the  sense  of 
viscous  asymmetry.  For  a  given  wavenumber,  easterly  shear 
dynamics  may  be  deduced  from  westerly  shear  dynamics  by 
interchanging  the  asymmetric  top  and  bottom  Ekman  layers. 

Mak  (1987)  in  a  two-layer  £-plant  model  found  that  a 
vertical  asymmetry  in  dissipation  qus//iat/ve{y  influences  the 
equilibration  of  nonlinear  quasi-geostrophic  unstable  baro¬ 
dinic  waves.  Especially,  he  found  that  the  sense  of  this 
asymmetry  is  a  crude/  factor  that  dictates  whether  the  equi¬ 
librated  state  is  generally  a  steady  multiple  wave  state  or  a 
vacillation.  It  is  uncertain  whether  Mak’s  results  arc  due  to 
the  sole  e/Tect  of  viscous  asymmetry  or  the  combined  effect 
of  that  asymmetry  with  0.  To  settle  this  question,  we  iso¬ 
late  the  effect  of  viscous  asymmetry  on  nonlinear  waves  from 
those  due  to  0  and  other  causes.  The  present  work  may  play 
the  role  of  a  bridge  between  our  previous  linear  analyses  and 
'.future,  more  realistic,  nonlinear  analyses. 
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2.  THE  MODEL 

We  use  the  model  of  Weng  (1990),  with  three  zonal 
wavenumbers,  each  having  two  meridional  modes,  in  wave 
field  and  four  modes  in  mean  flow  correction  (mfc).  It  is  an 
Eady-like  channel  model  having  asymmetric  Ekman  layers  at 
flat  top  and  bottom;  periodicity  applies  in  the  along  channel 
direction.  The  nondimensional  streamfunction  is  written  as 

d  *  -Ays  +  4>  +  d 

(basic)  (mfc)  (wave)  (1) 


The  nondimensional  governing  equation  and  boundary  con¬ 
ditions  are 


d,.  +  SV**»  0, 

i««0,  *t  y  -  0, 1,  _ 

hi  +  -  *h  -  2niv,4  +  J(hh)  -  Ahh)+ 

+  J(.hh)  +  Ahh)  =  0.  at  =  ±1/2. 


h.  +  *  °t 

d*=0,  at  y  w  0,1, 

ht  -  2s7<5„  +  4m)  ■  0. 


at  x  »  ±1/2, 


(3) 


where  \  )  stands  for  a  zonal  average  operator.  The  wave 
solution  is 


l  > 

?  =  53  H  + ec •  w 

tal 

where 

*ai  *  Ant  cosh(3/mx)  +  Bni  sinh(2/i,is)  (5) 

aad  e-e.  denotes  'complex  conjugate".  The  mean  flow  cor¬ 
rection  is 


w  ^[Af,  eoeh(2w,z)  +  Ar,sinh(2w,s)]sin(pry). 


The  nondimensional  parameters  and  the  notations  in  (l)  - 
(8),  as  well  as  those  need  in  this  paper,  are  listed  in  Table  1. 
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Abstract 

We  present  a  brief  overview  of  some  new  methodologies  for  making 
predictions  on  time-series  data.  These  ideas  stem  from  two  rapidly 
growing  fields:  nonlinear  dynamics  (chaos)  theory  and  parallel 
distributed  processing.  Examples  are  presented  that  show  the 
usefulness  of  such  methods  in  making  short-term  predictions.  It  is 
suggested  that  such  methodologies  are  capable  of  distinguishing 
between  chaos  and  noise.  Implications  of  these  ideas  and  methods 
in  the  study  of  weather  and  climate  are  discussed. 


1.  Introduction 

One  of  the  basic  tenets  of  science  is  making  predic¬ 
tions.  If  we  know  previous  behavior,  how  can  we 
predict  future  behavior?  The  approach  in  modern 
meteorology,  like  many  sciences,  requires  two  steps: 
construct  a  model  based  on  theoretical  considerations 
and  use  measured  data  as  initial  input.  Since  many  of 
the  underlying  theoretical  principles  in  meteorology 
are  known,  model  construction  has  been  and  contin¬ 
ues  to  be  a  primary  area  of  research  for  meteorolo¬ 
gists. 

Today’s  numerical  weather-prediction  models  for 
forecasting  tomorrow's  weather  (also  for  climate  pre¬ 
diction)  solve,  a  set  of  partial  differential  equations 
describing  fluid  flow  over  a  rotating  globe.  The  prob¬ 
lem  in  prediction  may  not  lie  here.  However,  as  was 
stated  by  Thompson  ( 1 957) ,  significant  problems  may 
arise  with  the  second  step,  where  measured  data  are 
used  as  initial  input  to  the  model.  Correct  specification 
of  initial  state  demands  the  measurements  of  vari¬ 
ables  in  a  three-dimensional  volume.  Routine  mea¬ 
surements  of  relevant  variables  are  taken  at  widely 
spaced  locations  providing  only  a  discrete  initial  state. 
The  spatially  continuous  differential  equations  simply 
cannot  operate  on  discrete  initial  input  (Farmer  and 
Sidorowich  1 987).  Because  of  this  inherent  forecast¬ 
ing  limitation  in  fluid-flow  problems,  we  are  motivated 
to  try  other  approaches. 

One  class  of  alternative  approaches  is  to  build 
models  directly  from  the  available  data.  For  these 
methods,  the  data,  given  as  a  time  series,  are  usually 
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considered  as  a  single  realization  of  a  continuous 
random  process  (see,  e.g.,  Pandit  and  Yu  1983).  As 
Farmer  and  Sidorowich  (1987)  point  out,  this  is  appro¬ 
priate  when  the  randomness  is  a  result  of  complex 
interactions  involving  many  independent  and  irreduc¬ 
ible  degrees  of  freedom.  Although  linear  methods  of 
analyzing  time  series  from  weather  and  climate  pro¬ 
cesses  have  had  some  success,  especially  in  regard 
to  relating  cause  and  effect  to  physical  phenomena, 
their  predictive  power  is  limited.  The  predictive  limita¬ 
tion  of  linear  methods  is  perhaps  related  to  their 
inability  to  model  feedback  dynamics  of  the  weather 
and  climate  systems  (Farmer  and  Sidorowich  1988: 
hereafter  referred  to  as  F388). 

In  the  last  decade,  advances  in  the  theory  of  dy¬ 
namical  systems  have  demonstrated  the  existence  of 
dissipative  systems  whose  trajectories  that  depict 
their  asymptotic  final  states  are  not  confined  in  limit 
cycles  (periodic  evolutions)  or  tori  (quasi-periodic  evo¬ 
lutions),  but  in  sets  of  the  total  available  phase  space, 
which  are  not  topological.  These  sets  are  fractal  sets 
and  are  often  called  strange  attractors.  The  corre¬ 
sponding  dynamical  systems  are  called  chaotic  sys¬ 
tems  and  their  trajectories  never  repeat.  Thus,  their 
evolution  is  aperiodic  but  completely  deterministic. 
Because  the  evolution  is  aperiodic,  any  signal  mea¬ 
sured  from  a  chaotic  dynamical  system  “looks’  quite 
irregular  and  exhibits  frequency  spectra  with  energy  at 
all  frequencies  (broadband  spectra)  similar  to  those  of 
random  signals  (see  Tsonis  and  Eisner  1989  for  a 
discussion  of  chaos  and  weather).  Another  important 
property  of  chaotic  dynamical  systems  and  theirstrange 
attractors  is  the  divergence  of  initially  nearby  trajecto¬ 
ries.  Due  to  the  action  of  the  attractor,  the  evolution  of 
the  system  from  two  (or  more)  nearby  initial  conditions 
will  soon  become  quite  different.  Since  the  measure¬ 
ment  of  any  initial  condition  is  subject  to  some  error, 
such  a  property  imposes  limits  on  long-term  predic¬ 
tion.  Nevertheless,  for  a  short  time,  nearby  trajectories 
may  not  diverge  significantly,  and  thus,  even  though 
each  evolution  might  be  quite  complex,  knowledge  of 
the  dynamics  and  especially  of  the  structure  of  the 
attractor  (e.g.,  dimensions.  Lyapunov  exponents)  may 
prove  beneficial  to  the  goal  of  short-term  predictions. 

For  a  system  containing  many  irreducible  degrees 
of  freedom,  the  linear  statistical  approach  is  probably 
as  good  as  any  and  may  even  be  optimal  (FS88).  if. 
however,  the  irregular  behavior  is  a  result  of  low- 
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dimensional  chaos,  nonlinear  models  ought  to  be  able 
to  do  much  better  at  prediction  than  simple  linear 
models.  In  fact,  since  chaos  does  npt  occur  unless  the 
system  is  to  some  extent  nonlinear,  nonlinear  models 
are  necessary  to  approximate  chaotic  dynamics. 

The  purpose  of  the  present  paper  is  to  outline  some 
recent  advances  in  modeling  time  series  and  to  dem¬ 
onstrate,  through  the  application  of  a  particular  tech¬ 
nique,  their  usefulness  in  making  short-term  predic¬ 
tions  over  standard  autoregressive  models.  The  paper 
is  not  intended  to  be  definitive;  rather,  it  serves  as  an 
interim  report  on  time-series  modeling  efforts  currently 
being  explored  in  the  physics  and  applied  mathemati¬ 
cal  communities.  No  attempt  is  made  to  sort  out  the 
particular  advantages  and  disadvantages  of  the  vari¬ 
ous  methods  mentioned.  More  details  concerning 
particular  methods  and  applications  are  given  in  some 
of  the  references  provided. 

The  problems  of  weather  and  climate  forecasting 
offer  a  unique  arena  for  testing  and  developing 
nonlinear  prediction  algorithms,  not  only  because 
current  numerical  weather-prediction  models  are  lim¬ 
ited  to  some  extent  in  their  prediction  capabilities,  but 
because  long-term  reliable  observational  records  have 
recently  been  made  available  for  climate  research. 
Diagnostic  studies  with  these  data  utilizing  nonlinear 
prediction  schemes  are  a  required  step  in  the  direction 
of  understanding  and  quantifying  the  complexity  of  the 
global  weather  and  climate  systems. 

2.  Nonlinear  prediction 

The  term  “nonlinear  prediction"  covers  a  broad  spec¬ 
trum  of  methodologies.  Our  focus  here  is  on  dynamical 
state-space  models.  The  two  components  of  such 
models,  determinism  and  state-space  representation, 
can  be  considered  separately.  In  fact,  more  familiar  in 
meteorology  are  two  variants  of  these  models  contain¬ 
ing  one  or  the  other  component.  The  analog  method, 
suggested  by  Lorenz  (1969),  while  not  strictly  limited 
to  time-series  data,  is  essentially  a  deterministic  non- 
state-space  model.  It  is  based  on  the  idea  of  finding  a 
historical  weather  pattern  (analog)  that  closely  re¬ 
sembles  the  current  weather.  The  evolution  of  the 
historical  analog  provides  a  model  for  the  evolution  of 
the  present  weather.  Although  the  method,  as  tried  by 
Lorenz,  was  successful  in  estimating  practical  limits 
on  atmospheric  predictability,  it  was  largely  unsuc¬ 
cessful  for  operational  forecasting  due  to  the  lack  of  an 
adequate  history  of  large-scale  weather  patterns.  With 
longer  data  archives  now  available  and  a  focus  on 
limited  areas,  the  analog  method  has  reemerged  and 
appears  to  hold  promise  for  weather  forecasting  (e.g., 
Van  den  Dooi  1989;  Toth  1989). 


Another  class  of  related  nonlinear  models  are  the 
threshold  autoregressive  (AR)  models  (Tong  and  Urn 
1980;  Tong  1983).  These  models  rely  on  a  state- 
space  representation  but  are  essentially  statistical, 
having  a  deterministic  component  limited  to  a  single 
variable.  While  this  is  a  considerable  improvement  (for 
modeling  chaos)  over  the  strictly  statistical  linear  rep¬ 
resentation  of  the  traditional  AR  models,  it  may  not 
provide  enough  nonlinearity  for  geophysical  signals  in 
general.  Recently,  Zwiers  and  von  Storch  (1 990)  have 
shown  that  such  models  are  quite  useful  for  modeling 
the  Southern  Oscillation. 

Building  a  dynamical  state-space  model  from  time- 
series  data  requires  two  steps:  finding  an  appropriate 
state-space  reconstruction,  and  then  cnoosing  a 
nonlinear  representation  that  maps  visited  regions 
into  regions  not  yet  visited  in  the  reconstructed  space. 
The  state  space  can  be  replaced  by  the  phase  space 
using  the  method  of  delays.  This  is  done  by  taking  a 
scalar  time  series  x(t)  and  its  successive  time  shifts  as 
coordinates  of  a  new  vector  time  series  given  by 

X(0  *  (x  (0,  x  ((+t) . x  (( +  (n  -1  )t)J,  (1 ) 

where  n  is  the  dimension  of  the  vector  X(t)  and  t  is  a 
time  delay  taken  to  be  some  suitable  multiple  of  the 
sampling  time  Af  (see  Packard  et  al.  1980;  Takens 
1981).  Thus,  for  an  n-dimensional  phase  space,  a 
“cloud"  of  points  wilt  be  generated.  From  this  cloud  the 
various  dimensions  and  exponents  can  be  calculated. 
The  proper  choice  of  x  to  obtain  a  suitable  reconstruc¬ 
tion  has  been  the  subject  of  considerable  debate.  In 
principle,  %  can  be  any  length.  However,  if  it  is  too 
small,  then,  in  general,  x(t)  will  be  nearly  equal  to  *(( 
+  i)  and  not  enough  separation  will  exist  between  the 
chosen  coordinates.  If  the  dynamics  take  place  on  an 
attractor  of  dimension  N,  then  it  is  necessary  for 
determinism  that  n  z  N  (i.e.,  the  attractor  must  be 
embedded  in  at  least  its  dimension,  otherwise  it  fills  the 
embedding  space,  thus  behaving  like  a  random  pro¬ 
cess)  (FS88).  For  proper  reconstructions,  Takens 
(1981)  showed  that  n  -  2Af  + 1  is  sufficient  at  least  in 
principle. 

There  are  other  ways  to  construct  a  phase  space. 
The  use  of  derivatives,  whereby  the  coordinates  of  the 
phase  space  are  successively  ordered  higher  deriva¬ 
tives  instead  of  discrete  time  shifts,  is  an  alternative.  In 
fact,  this  is  the  underlying  intuitive  concept  of  the 
method  of  delays,  but  it  is  not  recommended  in  prac¬ 
tice,  except  for  perhaps  extremely  dean  data,  since 
differentiation  amplifies  noise.  A  better  alternative  is 
the  method  suggested  by  Broomhead  and  King  (1986), 
whereby  the  Karhunen-Loeve  principal  value  decom¬ 
position  is  applied  to  the  vector  time  series  in  Eq.  (1 ). 
The  procedure,  called  singular-spectrum  analysis 
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1.  Introduction 

PLTMG  is  a  software  package  for  solving  a  general  second  order  parameter  involved  system  of  el¬ 
liptic  partial  differential  equations  on  a  two  dimensional  domain  with  Dirichlet  and/or  Neuman  boundary 
conditions.  It  was  originally  designed  by  Bank  and  his  collegues(1979b,  1982b,  1985b,  1988,  1990a)  as  a 
tentative  program  to  study  the  theoretical  and  practiced  aspects  of  the  multigrid  iterative  method,  adaptive 
grid  refinement  and  error  estimation  procedures,  and  their  interation.  The  adaptive  procedure  embodied 
in  PLTMG  is  based  on  the  ideas  of  Babuska  (1986a,  1986b,  1986c),  Rheinboldt,  and  their  co-worker  (1986, 
1987).  The  resulting  systems  of  nonlinear  equations  are  solved  by  a  combination  of  the  approximate  Newton 
iteration  (See  Bank  and  Rose  1981,  1982)  and  the  hierarchical  basis  multigrid  iteration  (See  Bank, Dupont 
and  Yserentant,  1988,  Yserentant,  1985,  1986a,  1986b,  1986c).  Jhe  mesh  refinement  algorithms  and  data 
structures  used  in  PLTMG  originally  followed  the  ideas  of  Bank,  Sherman  and  Weiser  (1983).  The  a  poste¬ 
riori  error  estimation  procedure  used  for  the  adaptive  mesh  refinement  was  developed  by  Bank  and  Weiser 
(1986,  1985,  1981).  The  pseudo-arelength  continuation  procedure  of  PLTMG  results  from  the  joint  work  of 
Bank  and  Chan  (1986).  The  sparse  Gaussian  elimination  procedure  and  global  stiff  matrix  storage  used  for 
the  resulting  linear  sparse  system  of  solutions  was  the  joint  work  of  Bank  and  Smith  (1987).  The  deflation 
technique  used  for  dealing  with  instability  caused  by  singular  matrix  was  developed  by  Chan  (1984). 

The  purpose  of  this  report  is  to  present  the  implementaion  of  the  software  package  PLTMG  Edition 
6.0  on  the  Florida  State  University  (FSU)  SCRI4D  Vax  and  to  present  some  examples  of  PLTMG  run  on  the 
SCRI4D  Vax  with  the  IRIS  4000  color  graphics  system.  The  main  algorithms  embodied  in  the  subroutine 
PLTMG  are  explained  in  detail  along  with  the  relevant  theoretical  background  and  major  results.  Data 
structures  used  in  PLTMG  are  also  illustrated. 
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Abstract 


A  quasigeostrophic,  two-layer,  midlatitude  jS-plane  channel  model  is  used  to  explore  the 
effects  of  steady,  zonally  asymmetric  potential  vorticity  (PV)  forcing  on  baroclinic  topographic 
instability.  It  is  shown  analytically  that  the  alteration  of  the  zonal  mean  flow  produced  by  the 
interaction  of  a  resonant  wave  with  the  PV  forcing  may  enhance,  suppress  or  catalyze  the 
topographic  instability  depending  on  the  phase  of  the  PV  forcing  with  respect  to  the  topography. 
For  PV  forcing  that  is  maximized  upstream  (downstream)  of  the  topography,  an  otherwise  stable 
(unstable)  flow  can  be  destabilized  (stabilized).  If  the  PV  forcing  is  considered  analogous  to 
diabatic  heating  in  a  continuous  atmosphere,  land-sea  heating  contrasts  can  be  expected  to  play 
an  important  role  in  the  genesis  and  nonlinear  evolution  of  topographically  induced  instabilities 
at  middle  latitudes. 

The  weakly  nonlinear  evolution  of  the  system  and  the  combined  effects  of  topographic 
and  zonally  asymmetric  PV  forcing  on  intraseasonal  oscillations  at  middle  latitudes  are  briefly 
discussed. 


,2 


Introduction 


The  circulation  of  the  atmosphere  is  subjected  to  a  vast  array  of  mechanical  and  thermal 
forcing  mechanisms  that  span  a  wide  range  of  space  and  time  scales.  For  the  large-scale  motions 
in  the  atmosphere,  the  predominant  mechanical  and  (external)  thermal  forcing  mechanisms  are, 
respectively,  topography  and  land-sea  heating  contrasts.  Such  forcing  mechanisms  have  long 
been  recognized  as  playing  an  important  role  in  a  variety  of  atmospheric  phenomena,  ranging 
from  the  preferred  regional  development  of  the  middle  latitude  cyclone  waves  to  the  global  heat 
and  momentum  balance. 

In  recent  years  several  studies  have  demonstrated  the  importance  of  large-scale 
topographic  forcing  in  generating  and  sustaining  low  frequency  motions  in  numerical  models 
(e.g.,  Legras  and  Ghil  1985;  Nathan  1988;  Nathan  1989;  Mak  1989;  O’Brien  and  Branscome 
1990;  Tribbia  and  Ghil  1990;  Howell  and  Nathan  1990;  Chou  and  Loesch  1990)  and  in  rotating 
annulus  experiments  (Li  et  al.  1986;  Pfeffer  et  al.  1989;  Bemardet  et  al.  1990).  Although  it  has 
been  recognized  since  the  work  of  Chamey  and  Eliassen  (1949)  that  the  standing  waves  at  middle 
latitudes  owe  their  existence  in  large  part  to  the  zonally  asymmetric  mechanical  forcing  produced 
by  the  earth’s  topography,  it  was  Chamey  and  DeVore  (1979)  who  established,  using  a 
barotropic  model,  that  topographic  forcing  could  yield  nonpropagating,  exponentially  amplifying 
solutions.  These  amplifying  solutions  have  since  been  referred  to  as  topographic  or  form  drag 
instabilities.  Chamey  and  Devore  also  demonstrated  that  at  finite  amplitude  these  instabilities 
may  lead  to  a  multiplicity  of  equilibrium  flow  configurations  for  a  given  set  of  model 
parameters.  One  of  the  equilibrium  flows  corresponds  to  a  high  amplitude  wave  state  with  a 
relatively  weak  zonal  current,  and  thus  has  been  associated  with  anomalous  flows  in  the 
atmosphere.  Subsequent  studies  also  have  demonstrated  the  existence  of  multiple  equilibria  -  and 


3 


more  generally  multiple  regimes  -  in  simplified  barotropic  and  baroclinic  models  of  the 
atmosphere  (e.g.,  Chamey  and  Straus  1981;  Pedlosky  1981;  Legras  and  Ghil  1985). 

However,  because  topographically  unstable  modes  and  multiple  equilibrium  states  require 
near  resonant  flow  conditions  for  their  existence,  a  condition  that  may  be  difficult  to  satisfy  for 
realistic  zonal  winds  in  the  atmosphere,  their  importance  in  the  atmosphere  has  been  questioned 
in  recent  years  (Tung  and  Rosenthal  1985;  Lindzen  1986).  Nevertheless,  there  exists  evidence 
to  suggest  that  topographic  instabilities  may  indeed  play  an  important  role  in  the  low  frequency 
variability  of  the  atmosphere.  For  example,  Legras  and  Ghil  (1985)  and  Tribbia  and  Ghil  (1990) 
have  shown  using  high  resolution  numerical  models  on  a  sphere  that  topographic  instability  may 
be  an  integral  part  of  the  overall  dynamics  governing  intraseasonal  oscillations  in  the  Northern 
Hemisphere  extratropics.  To  further  explain  these  oscillations,  Jin  and  Ghil  (1990)  used  a 
midlatitude,  /3-plane  channel  model  to  examine  the  resonant  response  of  equivalent  barotropic 
flow  to  topography.  They  showed  that  low  frequency,  finite  amplitude  oscillations  emerge  when 
wave-wave  interactions  are  considered  in  concert  with  the  topographic  form  drag.  Those 
oscillations  were  shown  to  depend  critically  on  the  meridional  profile  of  the  zonal  mean  flow  and 
the  existence  of  a  dipole  shaped  resonance. 

Because  the  importance  of  topographic  instability  and  multiple  equilibria  to  low  frequency 
variability  remains  largely  unresolved,  it  is  especially  important  that  we  determine  the  robustness 
of  topographic  instability  to  the  presence  of  other  internal  and  external  forcing  mechanisms. 
With  this  in  mind,  we  extend  Pedlosky  (1981)  by  examining  the  role  of  steady,  zonally 
asymmetric  potential  vorticity  (PV)  forcing  on  baroclinic  topographic  instability.  Such  PV 
forcing  can  be  thought  of  as  mimicking  the  thermal  forcing  that  arises  in  a  continuous 
atmosphere  due  to  land-sea  heating  contrasts.  Thus  several  questions  immediately  come  to  mind: 
Most  obviously,  can  zonally  asymmetric  PV  forcing  affect  the  topographic  instability 
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mechanism?  If  so,  how?  How  important  is  the  phase  between  the  PV  and  topographic  forcings 
to  the  overall  dynamics  of  the  system?  Are  finite  amplitude,  low  frequency  oscillations  possible? 

Model  and  Governing  Equations 

To  address  the  questions  raised  above,  the  conventional  two-layer,  quasigeostrophic  model 
is  used.  Briefly,  the  model  consists  of  two  layers  of  homogeneous,  incompressible  fluid  of 
slightly  different  densities.  The  system  is  gravitationally  stable,  rotating  with  angular  speed  Q 
about  the  vertical  axis,  and  confined  to  a  zonally  periodic,  midlatitude  channel  of  width  L.  The 
beta-plane  approximation  is  used  to  model  the  effects  of  the  earth’s  sphericity,  i.e.,  f  =  2ft  = 
fo+£'y\  where  f0  =  2flsin<£0  is  the  Coriolis  parameter  evaluated  at  a  central  latitude  <f>0,  y'  the 
dimensional  latitude  coordinate,  and  0'  the  dimensional,  northward  gradient  of  f.  The  quasi¬ 
geostrophic  potential  vorticity  equation  can  be  written  in  nondimensional  form  as  (Pedlosky 
1987) 


dO ■  0Ur . 

-g  *  P-^  ♦  ■/(*>  Q)  -  <HV«>  ns)  -  -<■/*,  ♦  Df  <i) 

where  Q  =  V20j4-(-l)iF(i/'1-^2)  is  the  quasigeostrophic  potential  vorticity  and  ^(x,y,t)  is  the 
geostrophic  streamfunction  field;  the  subscript  j  =  1  or  2  denotes  variables  in  the  upper  or  lower 
layers,  respectively;  tjb  is  the  bottom  topography,  F  is  the  internal  rotational  Froude  number,  0 
is  the  planetary  vorticity  factor  and  Dj  represents  the  imposed  potential  vorticity  (PV)  sources 
and  sinks.  In  (1)  V2  *  d2/dx2+SVdf  and  J(a',b')  =  {dz'/dx)(db'tdy)  -  (dz'ldy)(db'/dx).  . 

Two  forms  of  dissipation  are  considered  in  the  model:  Ekman  damping  on  the  upper  and 
lower  boundaries,  measured,  respectively,  by  the  parameters  r,  and  r2,  and  thickness  damping, 
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Abstract 


* 


Jin  and  Gbil  (1990)  demonstrate  that  for  topographically  resonant  flow,  low  frequency 
finite-amplitude  oscillations  may  arise  from  wave-wave  interactions  and  topographic  form-drag. 
Their  model  is  extended  to  include  a  zonally  asymmetric  vorticity  source,  which  is  shown  to 
interact  with  the  perturbation  field  to  produce  zonally  rectified  wave  fluxes  that  dramatically 
alter  the  Hopf  bifurcation  from  stationary  to  low  frequency  oscillations.  The  frequency, 
intensity,  and  general  character  of  these  oscillations  are  shown  to  depend  crucially  upon  the 
phasing  and  relative  strength  of  the  forcings. 

1.  Introduction 

Topographic  (i.e.,  form-drag)  instability  is  among  the  mechanisms  proposed  as  playing 
an  important  role  in  atmospheric  low  frequency  variability  (e.g.,  Frederiksen  1983;  Lcgras  and 
Ghil  1985;  Nathan  1988,  1989;  Howell  and  Nathan  1990;  Frederiksen  and  Bell  1987;  Jin  and 
Ghil  1990;  Tribbia  and  Ghil  1990;  O’Brien  and  Branscome  1991).  As  originally  shown  by 
Charney  and  DeVore  (1979),  alteration  of  the  zonal  mean  momentum  by  the  topographic 
form-drag  is  crucial  to  topographic  instability  and  the  subsequent  development  of  multiple  flow 
regimes  at  finite  amplitude.  Therefore,  other  external  forcing  mechanisms  that  interact  with  the 
perturbation  field  to  produce  changes  in  the  zonal  mean  momentum  will  likely  play  an  important 
role  in  topographically  induced  instabilities.  This  point  is  underscored  by  Nathan  (1992),  who 
analytically  demonstrated  within  the  context  of  the  two-layer  baroclinic  model  that  steady, 
zonally  varying  potential  vortcity  (PV)  sources  can  interact  with  the  wave  field  to  alter  the  zonal 
mean  flow,  resulting  in  dramatic  changes  in  the  linear  topographic  instability  properties  of  the 
system.  However,  at  finite-amplitude  only  steady  state  solutions  were  obtained;  the  absence  of 
perpetual,  low  frequency  oscillations  was  attributed  to  the  exclusion  of  wave-wave  interactions 
within  the  model. 

The  goal  of  this  note  is  to  demonstrate  the  impact  of  steady,  zonally  varying  vorticity 
forcing  on  topographic  instability  and  finite  amplitude ,  low  frequency  oscillations.  We  use  the 
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quasigeostrophic,  barotropic  channel  model  of  Jin  and  Ghil  (1990;  hereafter  JG),  who  showed 
that  for  topographically  resonant  flow,  wave-wave  and  wave-mean  flow  interactions  are 
essential  in  producing  Hopf  bifurcations  from  stationary  to  low  frequency  oscillations  on  the 
intraseasonal  time  scale.  Their  results  are  consistent  with  those  of  Tribbia  and  Ghil  (1990),  who 
numerically  integrated  a  fully  nonlinear,  global  barotropic  model.  In  both  models  the  forced 
wave  field  was  shown  to  occasionally  drift  back  and  forth  over  the  mountain,  a  feature  not  found 
in  observations  or  GCM  simulations.  JG  suggested  that  the  absence  of  zonally  varying  thermal 
forcing  may  explain  this  discrepancy.  Here  we  extend  JG  by  introducing  a  barotropic  analogue 
to  this  thermal  forcing  and  demonstrate  its  importance  in  affecting  the  phasing  and  strength  of 
the  low  frequency  oscillation  regime. 

2.  The  model  and  governing  equations 

The  model  consists  of  a  viscous,  barotropic  fluid  on  a  zonally  periodic,  midlatitude 
0-plane  channel,  bounded  above  by  a  flat,  rigid  boundary,  below  by  a  bottom  topography,  and 
laterally  by  sidewalls  at  y  =  0,  x.  For  a  steady,  meridionally  sheared  basic  flow  U(y)  the 
nondimensional  equations  governing  the  quasi-geostrophic  dynamics  of  a  superimposed 
disturbance,  Y(x,y,t)  and  mean  field,  tj/(y,t),  are: 


r~r  Jv 


-  ah 


v>  +  (0-Uyy) ^  +  J(V,  V»  +  J(y.  hg)  =  -  yv V-^fy)]-  U^+qB. 


dx 


dx 


(2.1a) 


d  92y  82y  a_ 

9t  9y2  ^  3y2  9y 


(2.1b) 


with  =  Vyt  =  0  at  y  =  0,  x.  In  these  equations  bg(x,y)  is  the  bottom  topography,  qB(x,y)  is  a 
steady,  zonally  varying  vorticity  source,  'Ffy)  is  a  perturbation  zonal  momentum  driving,  and  f 
measures  the  Ekman  pumping  strength  on  the  upper  and  lower  boundaries. 

Equation  (2.1b)  indicates  that  changes  in  the  mean  flow  vorticity  are  due  to  damping. 
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